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ABSTRACT 


More is now known about the behavior of the iono- 
spheric densities and temperatures as a function of solar 
activity, solar zenith angle, season, and geographical loca- 
tion. The currently available mid-latitude ionospheric and 
neutral atmospheric models are coupled with an advanced 
three dimensional ray tracing program to see what success 
would be obtained in predicting the wave propagation condi- 
tions and to study to what extent the use of theoretical 
ionospheric models is practical. The Penn State Mk 1 iono- 
spheric model, the Mitra-Rowe D region model, and the Groves' 
neutral atmospheric model are used throughout this work to 
best represent the real electron densities and collision 
frequencies. The Faraday rotation and differential Doppler 
velocities from satellites, the propagation modes for long 
distance high frequency propagation, the group delays for 
each mode, the ionospheric absorption, and the spatial loss 
are all predicted. 

The prediction of the oblique incidence soundings 
shows that the group delays measured are bracketed by the 
predictions for a reasonable range of solar activities. 

Comparison of the predictions of the model with the 
current techniques shows that the maximum percentage dif- 
ference between the maximum usable frequencies for the paths 
considered are in all cases tested within 5%. 
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The prediction of the Faraday rotation and differ- 
ential Doppler shift shows that the departure of the pre- 
dictions from the measurements are no larger than the range 
of variation encountered in a given month. 

Particular attention was paid to the interpretation 
of ionospheric absorption in terms of D region densities by 
clearly showing the percentage of the absorption that takes 
place at different altitudes. The absorption figures cal- 
culated by the model suggests that the method can serve as 
a useful technique in checking models of the lower iono- 
sphere . 


CHAPTER I 


INTRODUCTION 

1 . 1 General Statement of the Problem 

To understand ionospheric radio wave propagation it 
is necessary to know the expected ray paths. These are 
needed to allow quantitative estimates to be made of the 
expected performance of systems employing propagation through 
the ionosphere. It is desirable to be able to predict the 
propagation paths, the fading modes, field strengths, and 
the relationship of these to solar and magnetic conditions. 
This requires a model that provides information on the 
electron and neutral densities as a function of location, 
season, time, and solar conditions as well as a ray tracing 
program to enable propagation conditions to be calculated 
for complete paths under any desired conditions. 

It was desired to select from currently available 
mid-latitude ionospheric and neutral atmospheric models 
those that are best able to represent the real electron 
densities and collision frequencies. It was then desired to 
see how these could be coupled together in a consistent 
manner without discontinuities in boundary regions. It was 
then of interest to couple these with an advanced three 
dimensional ray tracing program to see what success would be 
obtained in predicting propagation conditions. 

For many years it has been common to justify aero- 
nomic programs on the basis of their relevance to radio wave 
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propagation. Theoretically based ionospheric models are 
now becoming available. It is of the utmost importance to 
study to what extent the use of theoretical models is prac- 
tical and what improvements are vital to their application. 

Chapter I introduces the problem and presents a 
comprehensive review of the relevant literature. 

Chapter II describes the method of calculation in 
predicting the wave propagation conditions by combining the 
theoretical ionospheric model with the ray tracing program. 
The coupling of the ionospheric and neutral atmospheric 
models is discussed in section 2.1 to obtain the continuous 
profiles of electron density and the collision frequency. 
Section 2.2 describes the techniques for developing and 
interpolating the ionospheric models and magnetic field 
model to supply the necessary parameters to the ray tracing 
program so that the ray path for any desired conditions can 
be calculated. Initial information needed to run the program 
is also described briefly in section 2.3. 

In the first section of Chapter III, predictions and 
observations for satellite propagation are compared under a 
variety of geographical, solar, and seasonal conditions. The 
horizontal gradient of the electron content along the satel- 
lite path is also discussed. In section 3.2 the prediction 
of the ray paths is described and the predictions of the 
oblique incidence soundings and the maximum usable frequency 
are compared with the experimental observations. In the last 
two sections the absorption of obliquely propagated high 
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frequency waves is re-examined and compared with the stan- 
dard absorption measurement technique. The spatial loss due 
to the power spreading is also predicted to see its contri- 
bution to the path loss. 

Chapter IV describes the prediction of the absorp- 
tion of the high frequency radio waves and its comparison 
with the experimental observations for a variety of condi- 
tions. The variation of absorption with solar zenith angle, 
ground range, and solar activity is described and an empiri- 
cal equation is derived. 

Chapter V discusses the conclusions of this study 
and compares the results with other studies. Conclusions 
given in the first five sections relate to each of the 
specific problems posed in section 1.3. Suggestions for 
future work are given in section 5.6. 

Appendix 1 describes the estimates of the initial 
azimuth angle and a great circle path for the location of 
the propagation path. Appendix 2 describes the main com- 
puter programs used in this study. 

1 . 2 Previous Related Studies 

1.2.1 Historical Summary of Ionospheric Prediction 

Because the ionosphere varies from hour to hour and 
from day to day, etc., it is necessary to have a knowledge 
of this variability in order to select the optimum frequency. 
The selection of the best sky-wave frequencies for communi- 
cation purposes over a given transmission path depends upon 
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the time of day, season of the year, phase of the sun-spot 
cycle, length of the path, and location of the path, as well 
as the conditions prevailing in the ionosphere over the path. 
The predictions are important for efficient frequency al- 
location, for the preparation of specifications for communi- 
cation equipment, and for the operation, planning, and design 
of high frequency communication systems. 

Current prediction techniques attempt to give values 
for several operational parameters. Those that have been 
commonly used are: 


MUF 

FOT 


M(3000)F2 


MUF (ZERO) F 2 
MUF (4000) F2 


The "maximum usable frequency" for a given 
circuit and a given mode of propagation. 

The "optimum working frequency" is defined 
as the value of the maximum usable frequency 
that is exceeded 90% of the days of the 
month. This frequency may be estimated for 
F2-layer propagation by multiplying the 
maximum usable frequency by 0.85. 

The maximum usable frequency for a 3000 km 
path divided by the F2-layer critical fre- 
quency for the ordinary mode. It is often 
called M(3000)F2 factor. 

The F2-layer maximum usable frequency for 
zero km. 

The maximum usable frequency of the F2- layer 
for a 4000 km path. 
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f Q F2 The F2-layer critical frequency for the 

ordinary mode . 

Predictions of the MUF are usually made for the - 
monthly median value, which is the value equalled or exceed- 
ed 50% of the days during the month at a specified time of 
day„ Many ionospheric observatories produce hourly values 
of the F2-MUF factor for reflections from the F2-layer. 

Many theoretical methods have been developed for the 
evaluation of the F2-MUF factor which is the ratio of the 
maximum usable frequency propagated to the critical fre- 
quency of the F2- layer for the ordinary ray* The parabolic 
layer theory, the non-parabolic layer theory, and the trans- 
mission curve technique are the most popular methods. 

The parabolic layer method depends upon the assump- 
tion that the F2-layer has a vertical electron density - 
distribution parabolic in shape. Appleton and Beynon (1940, 
1947) have developed, by assuming that the ionosphere has a 
parabolic electron density distribution, a direct method of 
calculating the maximum usable frequency , reflected at oblique 
incidence by a thick parabolic layer for both short-distance 
(plane earth) and long-distance (curved earth) transmission 
and explained the results graphically. They also showed how 
to determine the MUF and skip distance for various values of 
the half thickness y , the height of maximum electron density 

h , and the base height h , 
m o 

Newbern Smith (1939) gave a simple, rapid, graphical 
method for obtaining maximum usable frequencies and effective 
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reflection heights of radio waves, from vertical incidence 
measurements of the critical frequencies and the virtual 
heights of the various layers in the ionosphere . The method 
consists of the use of a "Transmission Curve", which is 
superimposed on the curve of frequency against virtual 
height, observed at vertical incidence. This transmission 
curve provided one of the most convenient and widely used 
methods for determining maximum usable frequencies. 

Another method is the so called non-parabolic layer 
theory (N(h) profile method) given by Vickers (1959). This 
method is based on empirical relations derived from opera- 
tional experience acquired over a long period of time. By 
producing a differential equation for the ray path which 
reduces, in the case of a parabolic electron density distri- 
bution, to that derived by Appleton and Beynon (1940) , and 
solving that equation numerically, oblique incidence propa- 
gation characteristics and MUF factors were deduced. These 
MUF factors were compared with those which would have been 
obtained had a parabolic electron density distribution been 
assumed or had a transmission curve been used. 

Up to this point most of the methods used in the 
determination of maximum usable frequencies for ionospheric 
propagation neglect the effects of the earth's magnetic 
field. Thomas, Haselgrove, and Robbins (1958a, b) and 
Thomas and Robbins (1958) have pointed out, however, that 
the inclusion of the magnetic field can appreciably affect 
the determination of the electron density distribution and 
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therefore the calculated MUF. Direct measurements of MUF 
(Wieder, 1955) have shown that the observed values are, on 
the average, about 3 percent higher than those obtained by 
the use of the transmission curves. Davies (1959) drew . 
attention to influence of the geomagnetic field on MUF 
calculations using a transmission-slider technique and gave 
suggestions for improving MUF calculations. 

As seen in the results of Hanson et al. (1958), 

Wieder (1955), and Silberstein (1958), the diurnal, seasonal, 
geographic variation of MUF factors and also a variation 
with sun-spot cycle could be seen. 

The first automated high frequency (HF) path 
prediction computer program was developed in 1957 for the 
U.S. Army Signal Corps. This was later modified in several 
versions. The first fully automated program in which the 
oblique transmission equations for parabolic layers were 
used was developed in 1966 (Lucas and Haydon, 1966) by ESSA's 
Institute for Telecommunication Sciences and Aeronomy which 
preceeded the Institute .for .Telecommunication Sciences (ITS) . 
In 1969 Barghausen et al. (1969) employed more extensive 
techniques which, though similar to the earlier ITS program, 
incorporated significant changes. 

One of the basic principles employed in, all iono- 
spheric prediction . is the relation of ionospheric character- 
istics to the sun-spot cycle. The actual prediction con- 
sists essentially. of first predicting the solar activity and 
then from the mass of data available, the corresponding 
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trends of seasonal, diurnal and geographical variations of 
the ionosphere characteristics. 

The basic ionospheric data used in the predictions 
are the E-, FI- and F2-layer ordinary wave critical fre- 
quencies (f 0 E, f FI, and f Q F2, respectively) and the 
F2-M3000 factor. The synoptic numerical coefficient repre- 
sentations of the ionospheric characteristics are fundamen- 
tal to all efficient HF computer prediction programs. These 
were first developed by ITSA and published in 1960 (Jones 
and Gallet, 1960) and developed further in 1962 (Jones and 
Gallet, 1962a; Jones and Gallet 1962b) and 1969 (Jones et 
al. , 1969). 

The basic form of the F2-layer predictions is a 
table of numerical coefficients defining a function which 
represents the world-wide and diurnal variations of iono- 
spheric characteristics like f Q F2, MUF(ZERO)F2, M(3000)F2 
median, and MUF(4000)F2 median. This function, referred to 
as the "numerical map" of the characteristics, has the form 
of a finite series of simple terms consisting of elementary 
functions of latitude, longitude, and time, each multiplied 
by the appropriate coefficient. A numerical map may then 
be written in the general form 

H 

fi(A,0,T) = a o (X,0)+ V [>j (A , 0) cosjT+bj (A,6)sinjT] 


where 
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k 

Ej 0,6) = ^a. Ck) F k (X,e) 
k=o 

k 

b.(x,e) = ^b. Ck) F k (x,e) 
k=o 

X = geographic latitude 
6 = geographic longitude 
T = universal time 

H is the number of harmonics for representing the diurnal 

variation. Here F k (X,0) denotes the orthonormal functions 

which are linear combinations of the geographic coordinate 

function G k (X,0). In most applications it is convenient to 

express fi(X,0,T) explicitly in terms of the G k (X,0). Then 

the new coefficient U , is used for defining the a.(X,0) 

s , x j 

and bj(X,0) in the form 

k 

Ej(X,6) - ^2i f k^k ^ ^ * j = ®»1»» • 

k=o 

k 

b.(X,0) = ^U 2j _ 1>k G k (X,0), j = 1,2,. . . ,H 

k=o 


Thus a numerical map fi(X,0,T) is defined by specifying the 
numerical coefficients to define the a^(X,6) and bj(X,8) 
in the above equations. 
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The "Ionospheric Predictions" issued monthly by the 
Ionospheric Propagation Predictions Group of ITS, ESSA 
Research Laboratories, Boulder, Colorado are based on this 
numerical map method. They are issued three months in 
advance and each issue provides the predicted sun-spot num- 
bers derived from a regression analysis based on the solar 
cycle, tables of numerical coefficients that define the 
functions describing the predicted worldwide distribution of 
f Q F2 and M(3000)F2, and maps for each even hour of universal 
time. of MUF(ZERO)F2 and MUF (4000) F2 . Davies (1965) and 
Ostrow (1966) explained the details of applying these 
numerical and graphical prediction maps and discussed the 
calculation of MUF and FOT for transmission by F2-layer 
alone, transmission by the regular layers (E-Fl, F2) , and 
transmission with the effects of the sporadic E E g for both 
distances less than 4000 km and distances greater than 
4000 km . 

Prediction of the long-term operational parameters 
of high frequency sky-wave telecommunication systems devel- 
oped by Barghausen et al. (1969) was an application of the 
numerical technique to the parabolic layer theory. The 
ionospheric parameters for the specific path such as critical 
frequency of the layer, height of the maximum electron den- 
sity of the layer, height of the bottom of the layer, and 
semithickness of the layer were obtained from the numerical 
world maps, and an empirical relationship between the maximum 
height of the F2-layer and M(3000) was used. Then the 
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parabolic layer theory was applied and the maximum usable 
frequency, take-off angle, and virtual height of reflection 
for- all frequencies were calculated'. This program predicts 
the long-term operational parameters such as maximum usable 
frequency (MUF) , and optimum working frequency (FOT) , based 
on the circuit reliability and service probability. A 
method. for predicting estimates , of the short term' opera- 
tional MUF based upon measurements of the iocal -magnetic 
activity in the vicinity of the circuit was also incorporat- 
ed into this program (Davies, 1970). 

There are still many areas where further study, 
modification, or extension of the prediction model and 
methods are desirable.- One may improve the prediction 
technique by using a ray tracing program which would be more 
adaptable when consideration must.be given to ionospheric 
'irregularities (tilts, spread F, sporadic E, etc.) that 
cause multipath, .mixed modes, . of fpath- scatter , asymmetrical 
paths, etc. Applying the ray tracing technique to a reason- 
able theoretical ionospheric model such as the Penn State 
Mk 1 Model should allow better predictions to be made. 

1.2.2 Propagation in the Ionosphere and Ray 
Tracing 

Since Kennelly and Heaviside (1902) independently 
suggested that long-range propagation of radio waves 'might 
be possible from England to America using the ionosphere as 
a reflector, numerous investigations of the ionosphere's 
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properties have been carried out for both practical and 
scientific purposes. 

The first quantitative study of radio- ionosphere 
interactions was made by Appleton (1932) who formulated the 
Appleton-Hartree equation for the complex refractive index 
and the polarization of plane waves propagating in a homo- 
geneous, anisotropic medium, such as the ionosphere. A very 
thorough discussion of them was given by Booker (1939, 1949), 
Ratcliffe (1959), Budden (1961), Kelso (1964), and Davis 
(1969) . 

The propagation of radio waves through the iono- 
sphere can be described quite simply if the index of 
refraction does not change appreciably over the distance of 
a wavelength. In this case the W.K.B. solutions which 
Wentzel, Kramers and Brillouin (1920's) used extensively in 
quantum mechanics, often produce very good approximations to 
the analytical solution of the wave equations governing 
radio wave propagation in the ionosphere. A W.K.B. solution 
is, therefore, often said to be a mathematical expression 
of "ray" theory. The W.K.B. solutions play a most important 
part in the theory of the propagation of radio waves, and 
comprehensive reviews of their use may be found in texts by 
Budden (1961), Ginzburg (1961), Kelso (1964), and Davies 
(1969) . This is the theory on which ordinary ray tracing 
programs are based. 

There are, of course, varying degrees of sophisti- 
cation in ray tracing. All ray tracing methods rest on the 
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same foundation which can be stated mathematically in a 
variety of equivalent ways. For example, Snell's law and 
Fermat's principle are the underlying principles on which 
all ray tracing techniques are based. An extension of 
Fermat's principle leads to a corresponding extension of 
Snell's law and could be used to derive all other equations 
used for ray tracing. However, a variety of mathematical 
approaches can be chosen to meet their own logical or 
computational purposes. Extensive descriptions of programs 
have been given by Croft (1969) , Kelso (1968) , and Jones 
(1966, 1970) . 

The ray tracing program used for this work is the 
Three Dimensional Ray Tracing Computer Program developed by 
Jones (1966) at the National Bureau of Standards. This ray 
tracing program has been in use now for over three years by 
people scattered all over the world. The method of tracing 
rays based on the principles of Hamiltonian Optics developed 
by Haselgrove (1955) is used in this program. The method 
involves the solution- of a set of six simultaneous differ- 
ential equations in terms of the three position coordinates 
and three variables indicating the direction of the wave 
normal . 

The program calculates ray paths by numerically 
integrating the following six differential equations similar 
to those described by Haselgrove (1955) . 
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|f- = (V r -Real{n|^-})/Real(nn’} 

v r 

|| = (V 0 -Real{n|~--})/(rReal{nn , }) 

0 

H = (V, -Real{n|^— })/(rsin0Real{nn' }) 

dV r Real{n|—} ,, 

3t~ = Rea'Knn 1 } * v e 3t + V in6 3t 


dv 0 1 I Real{n|f} 


■ v sl* *V°* e &. 


dV # j /ReaUnff) . ndr „ de 

3t ‘ rsTn? UeaHnnM ' V * slnS 3t ' r V° S 3t 


where the variables r, 0, and <j> are the spherical polar 
coordinates of a point on the ray path; V r , V Q , and are 
the components of the wave normal direction in r, 0, and <j> 
direction, normalized so that 


V r 2 + V e 2 + V^ 2 = Real(n 2 } 


n is the complex phase refractive index and n* is the complex 
group refractive index. These are expressed as follows. 


n = 1 


1 - jZ - 


4(l-X-jZ) 


2 !/2 
+ Y L } 


e mw 
o 
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N 


Y 

£ 

B 

c 

e 

Z 

v 

n' 


H 


electron density 
mass of an electron , 
free space permittivity 
2nf t angular wave frequency 
wave frequency 
plasma wave frequency 

longitudinal component of the geomagnetic field 
parameter Y along the wave normal direction 
transverse component of Y 

H B o 

mw 

1 I e I 

Tif — in — * electron gyro frequency 

the flux density of the earth's magnetic field 

charge of an electron 

v/w 

electron collision frequency 

n+f^, the complex group refractive index 


The group path length is used as an independent 
variable t in the above six simultaneous differential 
equations, whereas Haselgrove uses phase path length. The 
group path is defined by 


t 



n' cosads 


path 


where a is the angle between the ray direction and the wave 
normal direction. 
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Jones (1970) explained that the six simultaneous 
differential equations give the same path as that obtained 
by the equations of Haselgrove if there are no collisions. 

If there are electron collisions, so that n and n' 
are complex, then ray paths calculated using the equations 
of Haselgrove (1955) satisfy Fermat's principle, 

6 J' Real{n}cosads = 0 
path 

whereas, ray paths calculated using the above six equations 
satisfy 


6 J y Real(n 2 Jcosctds = 0 
path 

which is a better approximation to the ray path in a lossy 
medium. Jones (1970) mentioned that this ordinary iono- 
spheric ray tracing program works well at HF (3-30 MHz) . 

For the extension to lower frequencies, Jones also described 
the ray tracing in complex space to calculate the correct 
path in a lossy medium and compared the results with the 
phase integral method and the full wave solutions. 

In addition to the six basic equations necessary to 
calculate the ray path, the program also integrates 

dP _ Real(n 2 } 

Reallnn ' j 
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to give the phase path, and integrates 


dA _ 10 2 tt f | imagin^ } | 

3T ~ log a i0 C Realinn ' 1 

V 


to give the absorption in db. One of the computer outputs 
for the propagation from Annapolis to State College is 
shown in Table 1. 

The refractive index equations used in this ray 
tracing program are based on the Appleton-Hartree formulas. 
The following quantities are calculated in the program. 


n 


nn' 


3n 3n 
n 3P 


n 


3n _3n 


3 q : 


n 


3V r » 


n 


3n 

3 V 



and 


_3n 

n 37 


In order to calculate the above quantities, the following 
quantities are provided: 


X, 


3X 
3r ’ 


ax ax 

30* 3<f>* 


, 3X 
and ^ 


from the electron density model, 


Y 3Y 3Y 3Y Y ^r ^r ^r 
* 3r» 39’ 3 <J) * I r* 3r » 36 * 3<j> * 

3Y fl 3 Y fl 3 Y fl 3Y , 3Y. 3Y. 

Y e —> “• V TT’ 39^’ and 3^ 


from the magnetic field model, and 



Table 1: One of the Ray Tracing Computer Outputs for the Propagation 

Between Annapolis, and State College 

Ordinary 00A Theoretical Profile for D- and E-Region. Date: 08/16/7Z 
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7 3Z_ 3Z 3Z 

» 3r» 3 6 * W 

from the collision frequency model. 

Therefore, when using this program, one must specify 

* 

the ionospheric models wui^ define electron density, 
collision frequency, and the earth's magnetic field as a 
function of position in space. The models chosen for this 
work are described in the following four sections. 

For low frequencies (below about 1 MHz) and espe- 
cially for very low frequencies, the ionosphere cannot always 
be regarded as a slowly varying medium and full wave 
solutions are often required for many problems. Jones 
(1970) explained that even enough ray theory breaks down for 
low frequency (LF) radio waves, it can be best approximated 
to the full wave solution by extending Fermat's principle in 
complex space, which is the criteria of the ray path to have 
a minimum attenuation in addition to a minimum wave inter- 
ference, This is the basis for the complex ray tracing 
program. Jones (1970) stated that applying ray tracing in 
complex space to a plane wave incident on a plane stratified 
medium gives a result that agrees exactly with the results 
obtained by the phase integral method in the complex plane, 
and that agrees satisfactorily with the full wave solution 
above 30 KHz for all of his results shown. Below this 
frequency, the amplitude errors for vertical polarization 
and phase error for both vertical polarization and 
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horizontal polarization increase with decreasing frequency, 
and. therefore, the full wave solution is needed. 

In general solutions for full wave equations can 
only be obtained through numerical integration using 
computers, except for certain restricted ionospheric models 
and propagation conditions. There are usually three ways 
of approaching the full wave solutions based on the type of 
dependent variables for the integration. 

Clemmow and Heading (1954) and Budden and Clemmow 

(1957) obtained the coupled forms of the differential 

) 

equations governing radio wave propagation in the ionosphere. 

Budden (1955a, 1955b) described the methods for 
obtaining numerical solutions of the differential equations 
which govern the propagation of long and very long radio 
waves in the ionosphere at vertical and oblique incidence. 

In this method he obtained the first order simultaneous 
differential equations for plane wave propagation as derived 
from Maxwell's equations and the constitutive relations for 
the ionosphere and calculated the reflection coefficient 
matrix whose elements relate the total reflected wave to 
the total incident wave when the ionosphere is assumed to 
be horizontally stratified. The dependent variable is then 
the reflection coefficient matrix and not the wave field 
vector. 

Barron and Budden (1959) resorted to another depen- 
dent variable, however, termed the admittance matrix which 
is related to the reflection coefficient matrix. The wave 
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admittance approach is used here, although the wave fields 
are variables preferred by others such as Pitteway (1965) . 

Pitteway (1965) extended Budden's technique to 
include the calculation of the wave fields and the polar- 
ization of long radio wave propagation in the ionosphere. 

A constraint was added onto one of the two separate indepen- 
dent solutions required in constructing the wave fields 
which are the dependent variables for the integration. 

Piggot, Pitteway, and Thrane (1965), using Pitteway's method, 
presented some of the characteristics of the waves reflected 
from the proposed ionospheric models and compared these 
results with ground measurements, made on low and very low 
frequency waves. 

Seliga (1965) derived the generalized equations 
which include collision frequency dependence on the electron 
energy and outlined the numerical techniques for their 
solution similar to that of Pitteway (1965) . 

In the derivation of the Appleton-Hartree theory, 
the collision frequency of an electron was assumed to be- 
independent of the electron velocity. However, laboratory 
experiments made by Phelps and Pack (1959) have shown that 
the collision frequency of electrons is very nearly propor- 
tional to the square of their velocity. Under these circum- 
stances Sen and Wyller (1960) , Budden (1965) , and Deeks 
(1966) have deduced the generalized magneto- ionic theory to 
describe the propagation of a plane electromagnetic wave 
through a magneto-ionic medium in which the collision 
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frequency of monoenergetic electrons is proportional to the 
electron energy. It was shown that the classical Appleton- 
Hartree formulas can be best approximated by using the col- 
lision frequency called the effective collision frequency or 
equivalent collision frequency. 

In the present work attention has been confined to 
those cases for which ray tracing is adequate and no attempt 
has been made to use full wave propagation programs with the 
models . 


1.2.3 The Penn State Mk 1 Ionospheric Model 
The Penn State Mk 1 Ionospheric Model has been 
described by Nisbet (1970a, b) . This model provides a com- 
puter program that allows a model of the complete electron 
density profile in mid- lat itudes to be calculated from 
120 km to 1250 km as a function of latitude, longitude, day 
number, time, and solar activity. 

Table 2 shows one of the computer outputs. Electron 
densities are listed for a given longitude, latitude, time 
(both universal time and local time), day number, and solar 
activity expressed in 10.7 cm solar flux. It also gives 
the values of electron content, the peak electron density, 
the percentage probability of observing sporadic E with a 
critical frequency greater than 0.5 MHz, the average value 
of the peak electron density when sporadic E is observed, 
and the standard deviation of the log-normal distribution. 

The Penn State Mk 1 Ionospheric Model is constructed 


as follows: 
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E and FI Region 

1. An analytic fit to the CIRA (1965) neutral 
atmospheric model is used for all production and 
loss calculations. 

2. The FI region is treated using a photo- 
equilibrium model with 3 ionic constituents and 62 
ionizing radiation groups each of which varies with 
solar activity. Solar zenith angles are calculated 
based on altitude, location, and time. 

F2 Region 

3. The F2-layer uses theoretical models arranged to 
fit boundary conditions of the peak electron density, 
the height of the peak, and the electron density at 
1000 km. 

4. Values of the peak electron density are obtained 
using the C.C.I.R. (1967) Report 340 atlas of iono- 
spheric characteristics. 

5. Values of the height of the peak are obtained 
from a model based on reduced ionograms . 

Sporadic E Region 

6. The model is log-normal in electron density 

i n 

log— — 

P (N E > n) = 0 . 5P [1 -erf ( —) ] 

m s /To” 


where 

P = the probability of observing blanketing 
sporadic E 
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n = the average peak electron density when 
sporadic E is observed 

a = the standard deviation of distribution. 

The model assumes that the ion and neutral temper- 
atures are equal and calculates the electron temperature. 
Details of the construction of this model have been given 
by Nisbet (1970a) and a series of tables of electron 
densities for this model have been given by Nisbet (1970b) . 

To give an idea of the agreement of the model with 
observations, observations are shown on the left of Figure 1 
and predictions of the model are shown on the right. This 
is for a low solar activity. Figure 2 is another example. 
This time for a very high level of solar activity the model 
and observed ionosphere appear to behave very similarly, and 
the agreement in the FI region is of the order of 20%. 

1.2.4 The Mitra-Rowe Electron Density Model for 
the D Region 

For the D region, the model developed by Mitra and 
Rowe (1972) was used. This is a middle latitude model, and 
the magnetic latitude for this model lies between 30 and 
60 degrees. 

O 

It is a model which uses Lyman a, 1 - 8 A X rays, 

O 

1027 - 1118 A radiation, and galactic cosmic rays as the 
ionization sources. 

The 10c 7 cm solar flux is used to represent the low 
solar activity to moderate-high solar activity. Cosmic 
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rays, 1 - 8 A radiation, electron temperature, and the 
production of N0 + and excited state 0 2 + are scaled with the 
phase of solar activity. 

Six groups of ions are used to reproduce experimen- 
tally observed profiles of electron and ion densities. The 
ion chemistry of the D region was discussed in detail by 
Mitra and Rowe (1972). The simplified positive and negative 
ion chemistry used in the model is shown in Figure 3. This 
describes an ion chemical scheme consisting of a six-ion 
(four positive and two negative) model. 

For major neutral constituents up to 90 km atmo- 
spheric model by Groves (1971) was used. 

For minor neutrals, f^O, 0^, and 0 were obtained 
from Bowman et al . (1970), and 0 2 ( A g ) was obtained from 
Evans and Llewellgn (1970). C0 2 from Hays and Olivero (1970) 

and NO from Meira (1971) were also used in this model as 
were other minor neutral constituents. 

Figure 4 shows the electron density profiles for 
three different solar zenith angles for April 1971. The 
solar activity of this month corresponded to S^q y = 140. 

The model gives electron densities starting from 60 km. 

1.2.5 The Collision Frequency Model 

The neutral atmospheric model of Groves (1971) and 
the collision frequency calculations of Itikawa (1971) were 
combined to give the collision frequency profiles below 
120 km. The Groves Model gives atmospheric temperature. 


f 



Figure 3 Six-Ion Model for D Region 
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Figure 4 Electron Density Models. Below 120 Km 
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pressure, and density for the first of each month and for 
10° latitude intervals with values tabulated at 5 km 
intervals from 25 to 110 km. 

Using the equation given by Jacchia (1970), the 
number densities of N 2 , 0^ , and 0 were deduced as follows. 

n (N~) = .7811 

M o 

N (0) = 2Ap (I - ^-) 

M FT 

o 

n( 0? ) = Ap [i- (1+ . 21) - i] 

L M M 

o 

where 

- 23 - 1 

A = Avogadro's number = 6,02257 x 10 mole 

M = the sea-level mean molecular mass = 28,960 
o 

_ 3 

p = the density in g cm 

FT = the mean molecular mass at each altitude 
The variables p and M are also tabulated at 5-km intervals 
from 25 to 110 km in Groves' neutral atmospheric model. 

Once the number densities of major neutrals were 
obtained, the collision calculations of Itikawa were 
employed to give the following equation for the mono- 
energetic collision frequency of electrons with neutrals. 



32 


<v m > en = (1.292 x10~ 10 +2.135x10 12 T-5.216x 10" 17 T 2 ) n(0) 

+ (3.923x10" 10 +8.600x10" 12 T-8.428x 10' 16 T 2 ) n(N 2 ) 
+ (7.936xl0" 10 +4.292xl0" 12 T-2.858xl0" 16 T 2 ) n(0 2 ) 


where 

T = the absolute temperature of the electron. 

The mono-energetic collision frequency of electrons 
with neutrals calculated in this way for middle latitude is 
shown in Figure 5. These calculations are compared with 
experimental observations and other theoretical expressions. 

Figure 6 shows the seasonal and latitudinal 
variation of the collision frequency. At lower latitude the 
summer and winter profiles almost coincide. At higher 
latitudes the profiles for the two seasons differ appre- 
ciably . 

Above 120 km the Penn State Mk 1 Model was employed 
to give both the electron-ion collision frequency and the 
electron-neutral collision frequency. The electron colli- 
sion frequency with ions was obtained by assuming that 
electron, ion and neutral temperatures are equal. The 
collision frequency calculations of Itikawa were also used 
for the electron-ion collision frequency. 

In Figure 7 the mono-energetic electron collision 
frequency with ions and neutrals is shown. Because the 
electron density changes are much larger than the neutral 
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density changes, the ion collision frequency varies by a 
much higher percentage. The Penn State Mk 1 Model allows 
the electron-ion collision frequency to be calculated for 
each model used. 

1.2.6 The Earth's Magnetic Field Model 

The earth's magnetic field is of primary importance 
in the ionospheric propagation of radio waves. It brings 
about a variety of phenomena in the radio wave propagation. 
Anisotropy in the medium arises from the gyration of moving 
electrons and ions in a magnetic field. The frequency of 
the gyration is called the gyromagnetic frequency or simply 
gyro frequency f H . 

When the plane-polarized wave is incident upon the 
ionosphere, it splits into two waves which are propagated 
independently. One is called the "ordinary wave", and the 
other is called the "extraordinary wave". Under this 
circumstance the refractive index is more complicated than 
that calculated when the effects of the magnetic field are 
not being taken into account. The real part of the 
refractive index can be greater or less than unity; and 
waves of frequency less than the electron plasma frequency 
can propagate in the medium. 

For this work the magnetic field model described by 
Jensen and Cain (1962) was used. It is a series of coef- 
ficients for the spherical harmonic analysis of the magnetic 
field (Chapman and Bartels, 1962) as follows: 
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6 n 


X = 


X n tg n cos mQ + h n sin m61 Rn+2 


n=l m=0 


6 n 


Y = 


win r m • a 

n fg n sin mQ 


n=l m=0 


6 n 


h” cos me] R n * 2 


= 7 ^ Z n tg n cos me + h ” sin me) R °* 2 


n=l m=0 


where 


m 

n 

m 

n 

,m 


* 3 ? P n,m tcos 


mP Ceos d>) 
n,m 

sm <f> 


Z n = - (n+1) p n m (cos <j>) 


n,m 


R 


a 

a + h 


a = the radius of the earth 


h = height above sea level 

Here the function P (cos <t>) is a multiple of the associ- 

II | ill 

ated Legendre function defined as follows: 


P n,m tcos ' sln 


A r n-m, (n-m)Cn-m-l) n-m-2 x 

<t> [cos <p - — cos ♦ 


- (n r?^4)'Un } -ly(2n-l) (n ' m ' 5) 


n-m-4 


<p - . . . ] 
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X, Y, and Z are, respectively, the north, east, and 
vertical (down) components of the magnetic field vector F in 
Figure 8. The Greek letters 9 and <p are the longitude and 
the colatitude. 

The coefficient g™ and h™ used in this program are 
those given by Jensen and Cain (1962) . 

Having defined X, Y, and Z, the dip angle, the 
magnetic flux density B, and the gyro frequency f pj can be 
computed . 


/ 


- 1 

I = tan x 

lVx 2 +Y 2 i 

I 

in 

radians 

B = Yx^+Y^+Z 2 

B 

in 

- 4 

the unit of 10 wb/m 

f H = 2.8 B 


£ h 

in 

MHz 


One of the computer outputs is shown in Table 3 for an alti- 
tude of 300 km. This program provides the quantities I, B, 
and f H for any longitude, latitude, and altitude. 

1.2.7 Techniques for Measuring Ionospheric 
Abs orption 

There are several standard techniques for the mea- 
surement of ionospheric absorption of radio waves. Among 
these the Al, A2 , and A3 methods are principal techniques 
for measuring ionospheric absorption. 

Al Method This is also known as the pulse- 
reflection method and consists of the measurement of the 
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+ X (NORTH) 



Figure 8 Geomagnetic Convention for Describing 
the Elements of a Magnetic Vector: X, 

Y, Z, H, D, I, and F or B 



Table 3: One of the Computer Outputs Showing the Magnetic Field Variables 
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amplitudes of the pulse echoes reflected from the E region 
of the ionosphere. The virtual height for all measurements 
generally ranges from 95 km to around 115 km depending on 
the frequency used and the ionospheric conditions such as 
time, season, and solar activities. The absorption measured 
in decibels, L, is given by (Mitra, 1970). 

i A v 

L = - 20 log -S. - 20 log A (db) 
n n ' 

where 

1^ = daytime amplitude 

I = nighttime amplitude 
h^' = daytime virtual height 
h ' = nighttime virtual height 

An excellent summary of this technique has been given in the 
IQSY Instruction Manual No. 4 (Piggot et al . , 1963). 

As in the other techniques the question of separa- 
ting absorption contributions from the D- and E-layers has 
been a problem,, A brief discussion of this problem was 
given by Mitra (1970). 

Gnanalingam (1972) gave figures for the D- and E- 
region absorptions for the study of equatorial ionospheric 
absorption at Colombo, Ceylon. A substantial part of the 
absorption was found to take place near the reflection level. 
From his analytical method developed for the purpose of 
resolving absorption data, the percentage of absorption 
below 90 km and below 85 km was shown. 
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For both cases, the E region absorption was dominant. 
The percentage of absorption below 85 km was shown to be 45%, 
35%, 32%, 29%, 24% for the frequencies 1.33, 2.0, 2.2, 2.6, 
and 3.2 MHz, respectively. As the frequency increases the 
percentage of absorption below 85 km decreases. The varia- 
tion of absorption with solar zenith angle, season, and 
solar activity was also discussed. 

A2 Method This is also called the cosmic radio 
noise method or riometer method. This technique was first 
introduced by Mitra and Shain (1953) . The principle is 
that absorption is estimated by comparing the noise power 
received with the noise power that would have been received 
in the absence of the ionosphere at the same sidereal time. 
The absorption is estimated from the expression (Mitra, 

1970) . 

L = 10 log 10 (P o /P) db 

where 

P = the signal power received 

P Q = the signal power that would have been received 
in the absence of the ionosphere at the same 
sidereal time. 

Details can be found in the IQSY Instruction Manual 
No. 4. The choice of the frequency is such that the riometer 
is of use on frequencies above the F-layer critical frequency 
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only. A typical frequency, used in high latitudes to 
measure polar cap absorption, -is about 30 MHz. 

From the cosmic Radio Noise measurements Saha and 
Venkatachari (1970) discussed the results of the measure- 
ments of the collision frequency and the electron temper- 
ature. The temperature values were significantly lower than 
those obtained by either incoherent scatter techniques or by 
rocket probe methods. 

By calculating the absorption which a given riometer 
system would be expected to measure, Heisler and Hower (1970) 
discussed the interpretation of riometer measurements. By 
making a comparison between calculated and measured absorp- 
tion, they concluded that nighttime absorption reported by 
many authors was only an apparent absorption arising from 
errors in setting the level of the quiet-day curve. 

A3 Method This is the technique by which the 
absorption for the oblique radio wave path is commonly 
measured, usually at short distances. Recommended fre- 
quencies and distances for middle latitudes are 2-3 MHz 
and 100 - 400 km (Mitra, 1970), 

In this method a continuous -wave transmitter is used 
and the receiver output is registered on a pen recorder. 

The antenna is designed to select the ordinary component, 
since the absorption is smaller for the ordinary than for 
the extra- ordinary . 

Basically absorption by this technique is obtained 
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by comparing the daytime field strength with that of night- 
time, since nighttime observations are used for calibration. 

The absorption reduced to vertical incidence (J_) for 
frequency f^ (used at oblique incidence) < is 

given by (Schwentek, 1966) 

L J_ ( f A» °) = " cos 0 lL A (£ A» 

E d' l 'd 

= - cos a 20 log g~ 7 -p 

n n 

where 

a = the angle of incidence 
= daytime field strength 
E n = nighttime field strength 
= daytime half path length 
l ' n = nighttime half path length 
For the determination of the daytime ionospheric 
absorption, calibration was made depending on where the 
nighttime reflection occurs. Schwentek (1966) showed that 

is very close to unity and may usually be neglected, 
if calibration is possible using a sporadic E layer at night 
as a reflector. If the nighttime reflection is from the 
F region, calibration was made by multiplying the recorded 
field strength, E n (true nighttime field strength), by 
A'(F)/A’(E) and by a factor determined by the transmitting 
and receiving aerial polar diagrams for the direction in- 
volved. 



45 


Extensive studies with the A3 technique over a long 
period of time have been made in Germany (Schwentek, 1971) 
to determine the D region absorption. The method also has 
been used for the study of anomalous absorption in winter 
in middle latitudes. 

Besides the Al, A2, and A3 methods, there are other 
techniques called the £ m ^ n method and the sweep -frequency 
method. 

A critical examination of the techniques which have 
been used for measuring ionospheric absorption was made by 
Mitra (1970) in the special issue devoted to contributions 
from the Working Group on Electromagentic Probing of the 
Upper Atmosphere set up by the URSI . Comparison of differ- 
ent absorption techniques and their relative advantages and 
disadvantages were made in this issue and in the text of 
Davies (1969) . 

1 . 3 Specific Statement of the Problem 

The specific problems carried out in this work are 
outlined below. 

1.3.1. To develop a general theoretical approach to 
the radio wave propagation prediction techniques by com- 
bining the theoretical ionospheric model with the ray trac- 
ing program. 


1.3.2 To compare some oblique incidence ionospheric 
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soundings with the prediction of the propagation conditions 
for an ionospheric model. 

1.3.3 To compare the prediction of radio wave 
propagation using an ionospheric model with current pre- 
diction techniques for high frequency radio wave propagation 
for ground to ground paths. 

1.3.4 To compare the predicted satellite to 
ground propagation conditions with observations. 

1.3.5 To predict the path losses of high frequency 
radio waves and the contribution of absorption and spatial 
loss . 



CHAPTER II 


METHOD OF CALCULATION 


2 . 1 Coupling of the Ionospheric and Neutral Atmospheric 

Models 

In this work two ionospheric models are used. One 
is the Penn State Mk 1 Ionospheric Model representing the 
region above 120 km, the other is the Mitra-Rowe Model for 
the region below 120 km. Both models use different neutral 
atmospheric models. The Penn State Mk 1 Model uses an 
analytic fit to the CIRA (1965) neutral model; the Mitra- 
Rowe Model uses Groves’ seasonal and latitudinal atmospheric 
model and extends it to 120 km by extrapolation. 

In order to give the continuous gradient of electron 
density and collision frequency profiles at the boundary, 
the boundary conditions such as neutral constituents and 
temperature must be matched at 120 km. 

The determination of boundary values for neutral 
constituents and temperature has been a problem for many, 
years, and considerable controversy still exists (Mitra, 
1972) . 

The following attempt was made to make the boundary 
values for two models continuous in both electron densities 
and collision frequencies above and below 120 km. After 
obtaining the profiles for electron density and collision 
frequency from each model, a modification was made by hand 
to make the profile smooth at the boundary. For most cases 
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two models showed a slight difference due to the incom- 
patibility of the two International Reference Atmospheres. 

The parameters needed for the ray tracing are the 
phase refractive index, the wave normal direction, the 
derivative of position and wave normal direction with respect 
to group path, and the derivatives of refractive index with 
respect to position and the wave normal direction. These 
quantities are supplied from the neutral atmospheric model, 
the electron density model and the spherical harmonic 
expansion of the earth's magnetic field which is an integral 
part of the Penn State Mk 1 Ionospheric Model. Figure 9 
shows how the ionospheric model couples with the neutral 
atmospheric model to calculate the ray path using the ray 
tracing program. 

2 . 2 Technique for Developing and Interpolating the 

Ionospheric Models and Magnetic Field Models 
Once the propagation path has been located geograph- 
ically, as shown in Figure 39 in Appendix 1, profiles of 
electron density have to be developed along the path. In 
evaluating the ionospheric parameters, some of the pre- 
diction techniques use several control points (reflection 
areas) such as the midpoint of the path, E-layer control 
point, and F2-layer control point. Since the characteristics 
of the ionosphere at the point of reflection affect the MUF, 
that point is called a "control point". 



49 



Figure 9 Schematic Diagram Showing the Coupling of 
the Models With the Ray Tracing Program 
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When the ionospheric irregularities (such as tilts 
that cause off-path and asymmetrical path) exist, trans- 
mission along the great circle may not be possible depending 
on the direction of tilts. The geomagnetic field also can 
cause the deviation of the path from the great circle path 
depending on the direction of the field. 

Under these circumstances, it is necessary to obtain 
the electron density profiles all over the area involved in 
the transmission. 

Taking these effects into consideration, the grid of 
the electron density profiles between the transmitter and 
the receiver was made and shown in Appendix 1. 

The grid of the electron density profiles forms a 
three dimensional matrix, the element of which represents 
the electron density for a given longitude, latitude, and 
altitude. The interval between columns representing longi- 
tudes and the interval between rows representing latitudes 
can be chosen depending on the path length. 

This is easily done by obtaining the electron den- 
sity profiles for these longitudes and latitudes for a given 
day number, time (universal time), and solar activity 
expressed in 10.7 cm solar flux. 

Once electron density profiles at every intersection 
of the column and the row are obtained, the vertical and the 
horizontal gradient of the electron density at any point 
along the ray path can be evaluated simply by using linear 
interpolation as follows. The gradient of the normalized 
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electron density X in the direction of r, 6, and cf> are 


n . nfl + jU Jj n ,-,n( I J ,K) 
z z (1 + 1) - z(l) 


3n n(I . J ,K+1) - n(I , J ,K) 

3i = — 'xtm) ■ mi 




n = n(I,J,K) ♦ ff • ♦ fy • ^ ♦ |f • ** 


where 


Thus , 


. 3 

n(I,J,K) = the electron density in cm at Ith 

altitude, Jth longitude, and Kth lati 
tude 0 

z,X, and ip = altitude, latitude (north), and longi 
tude (east) . Latitude and longitude 
are both in radians. 


X = K • n/(f 2 ) 

3X y 3 n / f $ 2 ^ 

77 ’ K 3? /(f ’ 

3X v 3n /r r2s 

_ . . K • ^/(f ) 


3X 
3 <(> 


— / (f 2 ) 

alp ' ^ J 


where 
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X = normalized electron density 
f = wave frequency in MHz 
K = 80.5 x 10" 6 

r, 0, and <f> = spherical polar coordinates 
In a similar way, the grid of the earth’s magnetic 
field elements is made. 

From the magnetic field elements Z, X, and Y, as 
shown in Figure 8, the f^ r , f^ and f^ ^ which are r, 0, 

<j> components of the gyro-frequency f^, respectively are 
obtained with the proper signs by taking into consideration 
the use of spherical polar coordinates. Then the gradients 
of the normalized magnetic field vector Y (Y here is not the 
eastward magnetic field vector) and its components are 
obtained as follows: 


3Y 

r 

TT 


1 . 

r 


3Y 

i 

30 


1 f H . r d.J,K+l) - f H<r (I,J,K) 
7 ' 1 X (K+l ) - X(K) 


3 Y 1 f„ (I , J+l , K) - f r (I,J,K) 
3lT = 7 * * ijj(j+l) - vp (JJ 


f H r (I,J,K) 3Y 9Y \ 3 Y 

Y r = + TT- Az M'^J* AX + 7T * 


3Y 0 3 Y q 3Y q 


3 Y , 3 Y , 3 Y , 


Vii ± , V x 1 

3r~* T0~» IF"* 0* 9r * 3"0”* > x <j) 
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are obtained in the same way. Thus Y and its gradients are 
obtained as follows 


Y =-\/y 2 +Y fl 2 +Y, 2 
V r 0 <j> 


3Y 

3r 

i(V r 

3Y 

• r + 
3r 

K 

CD 

8Y S 

' Sr 4 * 

Y * 

3Y. 

3Y 

30 

7 tV r 

3Y 

r „ 
• TT + 

Y 9 

. !Ii. 

36 

Y 4> 

• 

30 J 

3Y 
3 <j> 


3 Y 

. —I + 

3<j> 

Y e 

. !Ii + 

30 

Y <J> 

9Y a 

. i) 

3(J> J 


These are the quantities which must be supplied to 
the ray tracing program. 

For the collision frequency model the gradient can 
be obtained in the same way as the gradient of the electron 
density was obtained. 

For the constant collision frequency and electron 
density profiles, the program provides the way to calculate 
the gradients using the cubic interpolation and extrapo- 
lation . 

2 . 3 Initial Information Needed to Run the Program 

Once the ionospheric model, the neutral atmospheric 
model, and the magnetic field model are supplied to the ray 
tracing program, one can start to run the program to cal- 
culate the ray path. First the frequency, the elevation 
angle, and the azimuth angle must be adjusted for the 
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desired transmission for a given path. To get the proper 
initial guesses to run the program, the following are 
considered. 

(1) The variation of the critical frequency along 
the great circle path. This will allow one 
to estimate the frequency. 

(2) The variation of the height of peak along the 
great circle path. This will enable one to 
see the possible asymmetry of the ray path and 
the reflection height. 

(3) The direction between the earth's magnetic 
field and the direction of the propagation to 
see the possible lateral deviation or asymmetry 
of the ray path. 

With these general ideas, the frequency can be 
approximately chosen using the "secant law". 


f 


ob 


f sec <j> 
v o 


where 

f k = frequency, of the wave incident obliquely on 
a flat layer 

f - the equivalent vertical frequency 
<j> 0 = the angle of incidence 

By incorporating the critical frequency of the layer, the 
height of the peak, and the given path length, a crude 
approximation can be made to calculate the angle of 
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incidence and the maximum usable frequency. If virtual 
height can be used instead of the height of the peak, more 
accurate values can be obtained. 

For long distance propagation, the earth's curvature 
must be taken into account. In this case Davies (1969) gave 
the following formulas. 

4[h Q +(D 2 /8a)] 2 +D 2 
4 [h Q + (D 2 /8a) ] 2 

sin(D/2a) 

1+ (h Q /a) -cos (D/ 2 a) 



where 

a = the earth's radius, 6370 km 

h Q = height of the reflection from a thin layer 

D = the distance (the ground range) 

By using the height of the peak of the F2-layer for h Q , <J> Q 

and f can be easily calculated, 
max * 

The ITS Ionospheric Prediction method was used to 

provide the initial estimates. For a given path less than 

4000 km, f and f can be found from the prediction map, 
c max r r 

and thus h Q and <j> Q can be approximated using the formula 
given by Davies. For a distance greater than 4000 km, path 
MUF and path FOT can be predicted based on the two control 
point method. 

With this information, the ray tracting program can 
be run using trial and error correction until the desired 
transmission paths are obtained. 



CHAPTER III 


RESULTS OF CALCULATION 


3.1 - Comparison of Satellite Propagation Measurements 

with Predictions 

One of the most basic propagation measurements is 
the electron content from observation of the Faraday rota- 
tion or the differential Doppler shifts for satellite to 
ground propagation. These measurements have been described 
by Bowhill (1958, 1962), Ross (1965), Ross et al . (1968), 
da Rosa and Garriot (1969), and Garriot et al. (1970). 

Such measurements have a very distinct relevance to 
the errors in navigational satellite systems, the errors in 
angular position of radio stars, and similar calculations 
when radio waves are propagated through the ionosphere at 
frequencies considerably in excess of the peak plasma fre- 
quency of the F2-layer. 

A knowledge of the electron content for a given path 
allows a correction to be made for such errors. As shown in 
Table 2 the Penn State Mk 1 Ionospheric Model also provides 
the value for the electron content for a given longitude, 
latitude, time, day number, and solar activity. In this 
section the prediction by model calculation is made and 
compared with the satellite propagation measurements. 

3.1.1 The Diurnal Variation of Electron Content 

Figures 10 and 11 show the prediction of the diurnal 
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Figure 10 


Diurnal Variation of Electron Content 
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Figure 11 Diurnal Variation of Electron Content 
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variation of the electron content at University Park, Penn- 
sylvania. Total electron content has been calculated for 
285 passes of the S-66 satellite observed at University Park, 
Pennsylvania (40.8°N, 78°W) between October 24, 1964 and 
March 16, 1965. A detailed analysis of these results has 
been given by Solomon (1965) . 

The dotted line in the figure represents the average 
of this measurement. The Faraday rotation technique was 
used to calculate the content. The solid line is obtained 
from the Penn State Mk 1 Model calculations. In February 
1964 and November 1964 the monthly mean value of 10.7 cm 
solar flux which was used for this calculation was approxi- 
mately equal to 73. The day to day variations in 10.7 cm 
solar flux were relatively small for this period. 

In Figure 10, as expected, both the prediction and 
the observation show that there seems to be a steady average 
increase in total electron content from about sun rise 
through the day to midafternoon. From sun rise through 
early afternoon ..in November 1964 the model calculation 
shows higher values than those obtained by satellite 
measurement as shown in Figure 11. 

Figure 12 shows another comparison with the measure- 
ments made at Huancayo, Peru (12.05°S, - 75.35°E), of the 
polarization rotation of the 54 MHz signals from the 
Satellite 1961 Omicron 1 from September 1961 to February 
1963 (Ross, 1966). This satellite was launched in a high- 
inclination orbit, approximately circular at about 1000 km 
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Figure 12 Average Diurnal Variation of Electron Content 
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altitude. The measurement of electron content, therefore, 
refers to the height range up to about 1000 km. The Penn 
State Mk 1 Model gives the values for the electron content 
up to 1250 km. Since the electron densities above 1000 km 
are very small, the electron content up to 1250 km is almost 
equal to the electron content up to 1000 km. 

The satellite measurements and the predictions by 
the model calculation show similar diurnal variations. The 
model calculations are about 301 lower than the satellite 
measurement at the diurnal maximum. Even though the model 
assumptions break down at this low latitude, the agreement 
for the electron content seems to be comparable to the day 
to day variations. 

Sometimes it is desired to know the diurnal variation 
of the electron content for different seasons with associated 
solar and magnetic activities. Figures 13, 14, and 15 show 
the diurnal variation of electron content at different sea- 
sons and different phases of the solar cycle at Stanford 
(37 . 5°N, - 122 °E) . 

The solid lines were obtained from the Faraday rota- 
tion measurement made with the geostationary satellites 
(Garriot et al. 1969). The cross marks were again obtained 
from the model calculations. The monthly mean values of 
10.7 cm solar flux in March, July, and December were used 
for spring equinox, summer solstice, and winter solstice, 
respectively. 

For spring equinox the agreement between the pre- 
dictions and the measurements is quite good except in 1968 




Figure 13 Diurnal Variation of Electron Content- at Stanford for Spring 
Equinox at Different Phases of the Solar Cycle 
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Figure 14 Diurnal Variation of Electron Content .at Stanford for Summer 
Solstice at Different Phases of the Solar Cycle 
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solar year. The prediction for 1968 solar year shows about 
10% lower values than the measurements during day time. 

For summer solstice considerable discrepancies of 
up to 40% between the predictions and the measurements can 
be seen through all four solar years. 

For winter solstice the prediction for 1967 solar 
year shows higher values than the measurement from the late 
morning through whole afternoon . 

3.1,2 The Day-to-Day Variation of Electron Content 

Figure 16 shows the. electron content observed for a 
satellite on north to south passes over State . College from 
October 1964 to September 1965. The contents correspond to 
different times of day as the satellite precesses. There is 
about a 30 to 40% difference between measurements and 
predictions from the middle of November to early December, 
whereas the predictions and the measurements are still good 
during May and June. 

The day-to-day variations are of the same order as 
the diurnal . variations of the electron content as shown in 
Figure 11 and Figure 14, The predictions are for monthly 
average conditions and it remains to be determined whether 
the departures are due to changes in the neutral atmosphere, 
the ionizing radiation or some other effects. 

3.1=3 Horizontal Gradient of the Electron Content 

Also of interest is the electron content gradient in 
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Figure 16 Day-to-Day Variation of Electron Content Through 
October 1964 to September 1965 
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the horizontal direction. In order to improve the satellite 
measurements, it is often necessary to predict what sort of 
horizontal gradient of the electron content is expected. 

Figures 17, 18, and 19 show three examples of the 
comparison of satellite observations and the models. It is 
apparent that even when good agreement is obtained between 
actual and measured values of the electron contents the 
horizontal gradients may depart greatly from the values 
given by models. This is hardly surprising because effects 
such as gravity wave perturbations can greatly change the 
horizontal gradient without major effects on the average 
total profile. Wind structures which are not always regular 
can also change from one night to the next. 

3 . 2 The Prediction of the Ray Paths, Oblique Incidence 
Soundings, and the MUF 
3.2.1 The Prediction of the Ray Paths 
An attempt was made to predict the possible ray 
paths between Washington (37°N, - 75°E) and London (52°N, 
o°E) for April 10, 1971 and solar activities similar to 
those of this day. 

Since there is a time difference in local time be- 
tween Washington and London, specific time was chosen to see 
the effect of the changes in the medium on the propagation. 

At 23:00 U.T. London will be in darkness and 
Washington in day light. This produces the electron 
density gradients shown in Figure 20. The ray paths are 
also shown for 2 hop-propagation and 3 hop-propagation and 
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Figure 18 -Daytime Gradient. of Electron. Content 
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Figure 19 Daytime Gradient of Electron Content 
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it can be seen that the distance between the paths is 
increasing as the reflection height increases. The reason 
for this is obvious from the map or plane view of the ray 
paths shown in Figure 21. The- peak electron density can be 
shown along the ray path and the maximum (peak) electron 
density decreases as shown in Figure 21. 

For the 2 hop-propagation, the ground reflection at 
2G is in the sea east of Newfoundland. For the 3 hop case 
the first ground reflection is almost exactly at Halifax, 
Nova Scotia and the second in the mid-atlantic south of 
Greenland. 

Figure 22 is another prediction of the propagation 
conditions - at 7:00 U.T. At London the ionosphere was 
already sunlit and Washington was in darkness. Here the 
electron density gradients are rather different from before. 

Figure 23 is a plane view of the same propagation 
conditions. As shown in the figure the. first ground hop 
3G has moved from Halifax, Nova Scotia out into the sea east 
of Newfoundland. The peak electron density gradients are 
almost at right angles to the ray path in this case. 

For these predictions, the- effect of the earth's 
magnetic field and the collision frequency of the electrons 
were neglected. 

3.2.2 The Prediction of the Oblique Incidence 

Sounding and the Maximum Usable Frequency 


Merely plotting ray paths is of limited interest and 
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Figure- 21 -Projection of -the Ray Path on the Ground 
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Figure 23 Projection of the Ray Path 


76 


an attempt was made to see what success might be obtained in 
predicting the shape of an oblique incidence ionogram. 

A virtual path- versus - f requency record is called an 
oblique . incidence ionogram and is obtained by photographing 
the B scan on slowly moving film or paper in which the prop- 
agation time delay between the time a pulsed sinusoidal 
signal is transmitted from a remote site and the time the 
signal is received is determined. 

Figure 24 shows an oblique incidence ionogram taken 
from NBS Report No. 31 Atlas of Oblique Incidence Ionograms 
along with the theoretical predictions. It is for a 1150 
km path from Sterling, Virginia (39.0°N, - 77.5°E) to St. 
Louis, Missouri (38.7°N, - 90.3°E) measured at about 9 
o'clock in the morning local time on April 24, 1952. Since 
the records of 10.7 cm solar flux for the time of the 
measurement were not available it was decided to bracket the 
estimated Zurich Sun-spot Number of 30 with 2800 MHz fluxes 
of 83 and 110 in the units of 10 ^ w/m^/Mz. 

The frequency scale is greatly expanded so the small 
differences are magnified. The corresponding take-off angle 
versus frequency curves for the predictions are shown in 
Figure 25- It can be seen that, as the wave frequency 
increases, the height of reflection of the low-angle ray 
increases very slowly, whereas that of the high-angle ray 
decreases more rapidly „ The frequency at which the low- 
angle ray and high-angle ray traces join is called the 




Figure 24 Oblique Incidence Sounding for Sterling - St Louis, 1150 Km 






Figure 25 Variation of Take-Off Angle with Frequency 
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junction frequency/’ and represents the maximum usable fre- 
quency (MUF) for the given path and condition. 

The MUF of the prediction for the solar activity, 

S^q 7 = 83 is 10.4 MHz and that for 7 = 110 is around 
12.7 MHz. This prediction shows therefore that the MUF for 
this path is proportional to the 1.4th power of the 10.7 cm 
solar flux. 

Another prediction was made for the path between 
Washington and London, and compared with those obtained by 
the ITS prediction method. The maximum usable frequencies 
for the distances 3000 km and 4000 km from the transmitter 
were predicted by the model and compared with those obtained 
by the ITS method. The results are shown in Table 4. 

Using the procedure based on the two control point 
method, the path MUF was also predicted and compared. It 
can be seen that the predictions for a reasonable range of 
solar activities bracket the observations in a satisfactory 
manner . 


3 . 3 Prediction of the Path Loss of High Frequency Radio 

Waves 

3.3.1 Theoretical Prediction of the Absorption 

Compared with A3 Absorption Measurements and 
Its Relation to the Ionospheric D- and E- 
Region 

The purpose of the work described in this section is 
to re-examine the absorption of obliquely propagated high 



Table 4; Comparison of MUF Model Calculation with 
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frequency waves and to determine in particular how A3 
absorption measurements are related to the D region of the 
ionosphere . 

Recent papers such as those of Schwentek. (1971) have 
presented extensive studies of D region absorption over the 
season and solar cycle based - on A3 measurements. 

First the prediction- was made for the path from 
Annapolis - (39°N, - 76.5°E) to State College (40.8°N, - 78°E) 
at 3.36 MHz. 

The attenuation figures were calculated using the 
three dimensional ray tracing program of Jones (1966) . 

For the earth's magnetic model, the constant gyro- 
frequency and the dip angle of the path midpoint were used. 
The gyro-frequency and the dip angle are 1.5 MHz and 71.3° 
respectively for the height range of interest. 

Since the attenuation figures of the measurement are 
not the absolute calculations of the daytime ionospheric 
absorption, but the daytime signal strength normalized to the 
nighttime- value , the predictions are also made to give the 
normalized path loss. Figure 26 shows the observed normal- 
ized path loss for this path compared with predictions based 
on the model. The current ionospheric models used here seem 
to give too little attenuation compared with the measure- 
ments, and the electron density at the lower heights is 
probably greater than given by the Mitra-Rowe model. For 
the path from Annapolis to State College, the absolute day- 
time absorption is greater than the normalized path loss by 



A3 ABSORPTION AT 3.36 MHz 
X 12 APRIL 1971 
o 16 APRIL 1971 
9 17 APRIL 1971 

-O-O- THEORETICAL PREDICTION, 
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Figure 26 Diurnal Variation of A3 Absorption with the Prediction for the 
Path from Annapolis to State College, at 3,36 MHz 
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about 5 db. The remaining figures in this section show the 
relation of the, daytime absorption to the ionospheric D- and 
E- regions. The differences between the absolute daytime 
absorption and the normalized path loss will.be discussed 
fully in the following section. 

Figure 27 shows the variation of absorption with 
height. It is apparent that the E region attenuation in 
this case is rather large compared with the attenuation 
lower down. It should also be noted that there is a con- 
siderable amount of attenuation close to the maximum altitude 
where the rays are traveling horizontally. 

Figure 28 shows the variation of the reflection 
height with time for the model. It is apparent that these 
are, as would be expected, quite large. As the solar zenith 
angle increases the reflection height also increases. 

Figure 29 shows the relative contributions from the 
region up to 85 km and also for the region above 85 km. It 
is apparent that the E region absorption dominates at all 
times in this instance. 

Another prediction for the A3 measurement from Ft. 
Collins (40.6°N, - 105. 1°E) to White Sands Missile Range 
(32.23°N, - 106. 5°E) is shown in Figure 30. The frequency 
used is 5 MHz and the ground range is 940 km. This figure 
shows different aspects from the previous figure. However, 
the attenuation above 85 km is still not negligible compared 
with D region absorption. 
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Figure- 29 Diurnal Variation o£ A3 Absorption in the D- and E-Regions 
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The model has been developed to allow rapid 
estimates to be made of the absorption of obliquely prop- 
agated radio waves in the high frequency band. It has been 
tested for a variety of frequencies, path lengths, and solar 
zenith- angles . Reasonable agreement with A3 measurements 
was also obtained. In many cases the E region absorption is 
not negligible and great care must be taken to interpret A3 
absorption measurements as indicating D region conditions. 

3.3.2 The Prediction of Loss Calculation Due to 
Power Spreading 

In high frequency propagation via the ionosphere, 
the path loss consists of several components, such as the 
spatial loss, the ionospheric absorption, and the focusing 
loss. For most cases the loss figures show that the spatial 
loss due to simple inverse square distance or dispersion of 
the energy is the main loss mechanism in radio wave prop- 
agation. In this section an attempt was made to see the 
contribution of the spatial loss to the path loss, and thus 
to deduce the absolute, value for the daytime ionospheric 
absorption by comparing the nighttime path loss with the 
daytime path loss. 

From the concept of the effective distance (which is 
the square root of the ratio of the actual defocused cross- 
secticn to the cross-section at unit distance from the 
source) , Davies (1969) combined the spatial loss and the 
focusing loss to give the following expression. 


where 


L t = the spatial loss combined with the focus- 
ing loss, and will be called the spatial 
loss hereafter. 

d = the effective distance, 
e 

A = the receiver cross-section. 

A u = the reference cross-section at unit 
distance (1 km) away from the source, 
a * the radius of the earth, 6370 km. 

0 = the elevation angle (take-off angle) . 

R and 4> = the ground range in km and in radians 
respectively as shown in Figure 31. 
dR/de = the derivative of ground range with 

respect to elevation angle and deter- 
mined from the ray tracing results in 
this work. 

An attempt was made to see what differences between 
daytime spatial loss and nighttime spatial loss might be 
obtained for both the Annapolis - State College path and the 
Ft. Collins - White Sands Missile Range path. 

For the nighttime model only the Penn State Mk 1 
Model is used to give the nighttime values. for the electron 
densities and the collision frequency. 



90 




Figure 31 Divergence of a Narrow Radio Beam 
Reflected from a Curved Ionosphere 
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From the results shown in Table 5, there are 5 to 
6 db differences between nighttime spatial loss and that of 
daytime for Annapolis to State College path while. the path 
from Ft. Collins to White Sands Missile Range shows only a 
2 db difference. 

Since the attenuation figures shown in Figure 26 
represent the normalized path loss, the differences of 
daytime and nighttime intensities have to be increased in 
this case by about 5 db to give a reasonable estimate of the 
daytime absorption assuming the antenna gain remains the 
same at the different zenith angles. 

Here the ionospheric absorption by nighttime E 
layer below 120 km is neglected due to the assumed model. 

When the nighttime E layer electron densities are large, 
there is some absorption depending on the path geometry 
(Wakai, 1972) and therefore care must be taken. 

To give some idea about the attenuation figures due 
to the spatial loss, an attempt was made to calculate the 
spatial losses for the same conditions on which the values 
of the- ionospheric absorption were made as will be shown in 
Chapter IV. Due to the complexity of the ray paths the 
values shown in Table 6 are somewhat scattered. 

3.3.3 Discussion of the Absorption Calculation 
Near the Reflection Level 

The propagation of radio waves through the ionosphere 
can be described quite simply if the index of refraction does 


Table 5: Calculations of Path Loss for Annapolis to State College Path and for 
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Table 6: Calculations of Spatial Loss 
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not change appreciably over the distance of a wave length. 
In this case the reflection and transmission coefficients, 
the height of reflection, and so on, can be determined by 
the W.K.B. Solution. 

The W.K.B. solutions at oblique incidence, when the 
effect of the earth's magnetic field is not included, are 
discussed in Chapter 9 of Budden (1961) . 

For horizontal polarization the W.K,.B. solutions 
are invalid if 




1 1 
7 T 

q 



>> l 


and for vertical polarization 






1 

2q 


4 

dz z 


1 

T 

q 




dz‘ 




>> 1 


where 

k = the free space propagation constant 
q = the normal component of the phase refractive 
index 

Thus the W.K.B, solutions are invalid for 1) long 

wave length (low frequencies), 2) near reflection levels or 

coupling levels, 3) electron density not slowly varying. 

The W.K.B. solutions also fail if the derivatives 
2 2 

dq/dz, d q/dz are large even when q is large. This would 
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be the case of reflection at a sharp boundary between two 
homogeneous media where the derivatives are infinite. Note 
that when the effect of collision is present, q is complex 
and never takes on the value of 2 ero. However, there is a 
region where |qj is small and the W.K.B. solution is not 
valid. The W.K.B. solutions also break down for the case 
of partial reflection and partial transmission as discussed 
in chapter 17 of Ratcliff (1959). 

The generalization of the W.K.B. solution can be 
found in chapter 18 of Budden. (1961) and in the paper by 
Budden and Clemmow (1957) , 

When the W.K.B. solutions are invalid, the phase 
integral method or full wave solution must be used. 

The purpose of this section is to discuss the 
absorption claculations by the ordinary ray theory and the 
full wave solution for the high frequencies of interest. 

For some N(z) or X(z) profiles, it is possible to 
make a change of variable in such a way as to obtain a 
standard form of differential equation such as Bessels 
equation or the Hypergeometric equation. After obtaining 
the set of independent solutions of the equations and 
applying boundary conditions to obtain the required solution, 
the reflection coefficient is determined by separating the 
solution into the upgoing and downgoing parts and taking 
their ratio. 


In this way Budden (1961) obtained the reflection 
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coefficient for the purely linear profile and the constant 
collision frequency as follows: 


R = i exp{ - 2 ikCh Q - | iC 3 k(l-iZ)/a} 


where 

C = cos 0j 

0j = the angle of incidence 
k = the free space propagation constant 
i = 

for the linear electron density profile, 

X = a(z-li) 
o 

and for the constant collision frequency profile, Z. Thus 
the absorption per one hop can be obtained as follows: 

L = - 8 , 68&n | R | (db) 

where 

| R | = exp{ - i C 3 kZ/a} 

For the purpose of comparison the following linear 
electron density profile is used to represent the practical 
profile for the region between 85 km and around 120 km. 

N = 1.314 x 10 3 ’(z- 85), electrons cm" 3 
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where 

z = the altitude in km 

Table 7 shows the absorption calculated from the 
ordinary ray tracing program used and the values from the 
full wave solution. 

As shown in the table the agreement between the two 
methods is not bad when the assumed value of the collision 
frequency is close to the collision frequency at the level 
where the reflection occurs. Since absorption near the 
reflection level is the main interest here, it will be 
wise to choose the value for the assumed constant collision 
frequency such that this value is close to the value near 
the reflection level occurring when this assumed collision 
frequency is included. 

Since the ray tracing program tries to calculate the 
absorption near the reflection level where |q| is very 
small, it tends to overestimate the absorption in this 
region. This is why the ray tracing program gives higher 
absorption than those obtained from the full wave solution 
as shown in Table 7. 

For waves of higher frequencies which are reflected 
at higher levels where the collision frequency is smaller, 
there should be smaller discrepancies between the two 
methods for the absorption near the reflection level. It 
should be remembered that the ray theory discussed here is 
the ordinary ray theory, not a ray theory applied in the 
complex space. Jones (1970) explained the details of 



Table 7: Comparison of Absorption Calculation by the Full Wave 
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complex ray tracing as discussed in section 1.2.2 of 
Chapter I. 

Recently Gnanalingam (1972) calculated 3.2 db, 1.5 
db, 1.3 db, 1.0 db, and 0.8 db for the phase integral 
corrections for the absorption of vertically incident 
waves of 1.33 MHz, 2.0 MHz, 2.2 MHz, 2.6 MHz, and 3.2 MHz, 
respectively. These absorption measurements were made at 
Colombo which is just outside the region where the quasi- 
transverse approximation is applicable to vertical 
propagation of the frequencies employed. (Dip angle 5°S, 
f^ = 0.093 Mhz at 100 km). 

It can be seen that there are lower values for the 
phase integral correction for waves of higher frequency. 

For the oblique incidence this phase integral correction can 
be expected to be less than that for vertical incidence of 
the same frequency. 

Recent papers, such as that of Jones and Foley 
(1972) , consider the reflection and coupling points to be 
located away from the real axis in the complex X-plane 
when the effects of electron collisions are included. The 
conditions for reflection and coupling were derived using 
the Booker quartic equation. The behavior of the reflection 
levels in the complex X-plane was also shown as the angle of 
incidence is varied. 

This is equivalent to a complex reflection height as 
discussed in Chapter 16 of Budden (1961) . When the phase 
integral is calculated in the complex z-plane, the phase 
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integral along the imaginary axis gives the additional con- 
tribution to the phase shift and the reduction of amplitude. 
These were not predicted from simple ray theory (Budden, 
1961). Thus the limiting factor of ordinary ray tracing 
programs is their failure to include the complex contribution 
to the phase integral (Jones, 1970). 

It would be worthwhile to compare the result of 
ordinary ray tracing and full wave solution for more practi- 
cal profiles of the electron density and the collision fre- 


quency. 



CHAPTER IV 


THE VARIATION OF ABSORPTION WITH SOLAR 
ZENITH ANGLE, GROUND RANGE, AND SOLAR ACTIVITY 

I 

4 . 1 The Diurnal Variation of Absorption 

The absorption of radio waves is a complicated 
function of frequency, geographic location, season, ground 
range, solar zenith angle, and solar activity. The inter- 
relationship between these controlling factors makes it very 
difficult to calculate the absorption of HF radio waves in 
a simple manner. 

As shown in Figure 26 the absorption varies with the 
solar zenith angle x- This diurnal variation is often 
expressed in terms of the power of cos x* In Figure 32 the 
solar zenith angle dependence of the predicted absorption 
for the path from Annapolis to State College is shown. In 
the prediction the absorption both in the morning and in 
the afternoon were restricted by the model assumptions to 
be the same for the same solar zenith angle. 

The predicted absorption can be approximated as 

follows 


L (prediction) = 32.1 cos''”'^ x db 

Many experimental results show a wide variation in 
the value for n, the power of cos x> ranging from 0.75 to 


ABSORPTION (db) 
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Figure 32 Variation of Absorption with Solar Zenith Angle 
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1.5 depending on the other controlling factors such as solar 
activity, season, etc. 

Figure 33 shows the variation of absorption with 
solar activity for the same path. This prediction shows 
that the absorption increases with solar activity expressed 
here in 10.7 cm solar flux. Assuming a linear dependence 
of absorption on the solar activity, the figure on the 
right hand side shows the variation of the absorption with 
10.7 cm solar flux for three different solar zenith angles. 
For the low solar activity, the dependence of absorption on 
cos x is similar to that for high solar activity shown in 
Figure 32. 

4 . 2 The Variation of Absorption with Ground Range 

An attempt was made to predict the absorption of the 
waves of higher frequencies in middle latitudes for spring 
equinox. For this purpose the path between Washington and 
London was chosen. Electron density and collision fre- 
quency profiles are assumed to be the same along the path. 
The magnetic field model of Jensen and Caine (1962) was used 
to give the earth's magnetic field elements. 

Figure 34 shows the variation of absorption with the 
angle of incidence. The predictions are compared with a 
formula given by Rawer (1952) shown below 

L a = 430(1 + 0.0035 R) cos 3/ 4 X sec<|> D (f ±f ) 
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Figure 34 Variation of Absorption with 
the Angle of Incidence 
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where 

L = the absorption in decibels per hop 

cl 

R = the sunspot number 
X = the solar zenith angle 
dp = the angle of incidence 
f = the wave frequency in MHz 

fp = the longitudinal component of the gyro frequency 
This formula represents the absorption in middle 
latitudes and the numerical values in this formula are 
based on analysis of vertical- incidence absorption measure- 
ments as Slough, England. 

As shown in the figure the empirical formula given 
by Rawer has a linear variation of absorption with sec dp, 
while the prediction shows the relation to be non-linear. 

This is because of the importance of the deviative absorp- 
tion in the calculations. 

For the same conditions the variation of the absorp- 
tion with ground range is shown in Figure 35 with the values 
of elevation angles indicated beside the dots. 

The figure shows the possible existence of high- 
angle rays and low-angle rays for a given ground range. For 
example three possible rays can be found to have the same 
ground range around 950 km. One with an elevation angle of 
16° is reflected at 127 km, the others with elevation angles 
of 20° and 30° are reflected at 137 km and 194 km respec- 
tively. These multiple modes are often seen in oblique 


lonograms . 
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Figure 35 Variation of Absorption with- Ground Range 
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For high elevation angles, above about 45° it can 
be seen that the absorption increases with the elevation 
angle. The reason for this is that the electron density 
increases as the wave approaches the height of the peak, 
while the collision frequency (mainly with ions) remains 
almost constant in this instance. 

For an elevation angle of 60°, the variation of 
absorption with height is shown in Figure 36. A large 
amount of deviative absorption can be seen at high altitudes 
while the height variation of absorption for low elevation 
angles, such as 20°, shows different relative contributions 
as shown in Figure 37. For both cases the total absorption 
per hop is about 9 db . 

Figure 38 shows the oblique ionogram for a path 
about 1000 km long and the values of absorption correspond- 
ing to low-and high-angle rays. In this case the absorption 
of low angle rays is higher than that of high angle rays. 
However, the variation in the magnitude of the absorption 
of low-angle rays is somewhat different from that of high- 
angle rays. Even though the absorption of low-angle rays is 
higher than than of high-angle rays for the same ground 
range, the spatial loss due to the energy spreading will be 
expected to be higher for high-angle rays than for low-angle 
rays. Details of the spatial loss were discussed in the 
previous chapter. 
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Figure 37 Variation of Absorption with. Height 
for an Elevation Angle of 20° 
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4. 3 The Variation of Absorption with Solar Zenith Angle 
and Ground Range for Two Different Solar Activities 
It was decided to attempt to derive an empirical 
equation relating the ground range, solar zenith angle, and 
10.7 cm solar flux. Calculations were made of the absorp- 
tion for three different solar zenith angles and a set of 
frequencies ranging from 5 MHz to 30 MHz. Two solar activ- 
ities S^q y * 70 and j = 140, which correspond approx- 
imately to the sun-spot numbers R = 0 and R = 100 respec- 
tively, were used to represent two phases of solar 
activities. The results are tabulated in Table 8. 

It should be remembered that within a small range 
of elevation angle there can be possible low-and high-angle 
rays even though the results are tabulated with a large 
interval of elevation angle. 

Even though the correct value for the absorption 
can only be obtained from the ray tracing, it is worthwhile 
to have- an empirical equation to give a quick estimate of 
the absorption. 

For this purpose the data obtained from the ray 
tracing were fitted to give the following equation. 

L = 2 00cos° ° 71 x(l + 0.0031 • S 1Q Q+0 . 0O41D) 

• {(f+f H ■ sin(I)) 2 +Cv o /2Tr) 2 }' 1 db 


where 


Table 8: Calculations of Ionospheric Absorption 
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L = total absorption per hop in db 
X = solar zenith angle 

Sio 7 = solar activity expressed in 10.7 cm solar flux 
D = ground range in km 
f = wave frequency in MHz 

= gyro frequency in MHz 
I = dip angle (at the midpoint of the path) 
v o /2tt = collision frequency at 90 km in MHz 
A nonlinear least square fitting method was used in this 
analysis. The standard deviation of the predicted values 
by this empirical equation from the data is around 3 db. 



CHAPTER V 


SUMMARY AND CONCLUSION 


5 . 1 Development of a General Theoretical Approach to 
Radio Wave Propagation Prediction by Combining a 
Theoretical Ionospheric Model with a Ray Tracing 
Program 

It was not the purpose of this work to develop new 
ionospheric models or ray tracing programs. Modification 
was made to allow the D region model to be matched to the 
bottom of the Penn State Mk 1 Model and to ensure that a 
collision frequency model was developed that was consistent 
with both the most up-to-date atmospheric model in the lower 
ionosphere and the Penn State Mk 1 electron densities in the 
F region. 

In contrast to current prediction techniques the 
method described in this report allows many properties of 
high frequency radio wave propagation to be predicted in a 
consistent manner. It has been shown that Faraday rotation 
and differential Doppler velocities from satellites, prop- 
agation modes for long distance high frequency propagation 
including the effect of longitudinal and transverse gradi- 
ents, group delays and azimuth and zenith angles for each 
mode, ionospheric absorption, and spatial loss including the 
effects of focusing can all be predicted. 

While the method requires large amounts of computer 
time for the ray tracing program, it does not require a 
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great deal of knowledge about ionospheric conditions or 
wave propagation on the part of an eventual user. The 
complexities in these areas are built into the program. The 
restrictions on the use of the technique depend on those 
of the ionospheric models and on the ray tracing program. 

The current level of knowledge of the ionosphere only 
allows average monthly ionospheric conditions to be pre- 
dicted. It is probably a fair summary of the work reported 
here to say that in each case studied, the departures of the 
predictions from the measurements were no larger than the 
range of variation encountered in a given month,. While the 
measurements of the densities and contents were within these 
limits, the estimation of gradients in the density was much 
more difficult. Large departures from average gradients 
were observed and would be expected, so that prediction of 
such things as the departure of the ray path from the great 
circle path and focusing should be regarded with great 
caution, A model, to account for these effects fully, will 
presumably require frequent updating with real time 
ionospheric data, 

5 . 2 Comparison of Oblique Incidence Ionospheric 

Sounding with the Predictions of the Propagation 
Condition for an Ionospheric Model 
An oblique incidence sounding for the 1150 km path 
from Sterling, Virginia to St, Louis, Missouri was compared. 
It was shown that the group delay and the relation of group 
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delay with frequency for both low- and high-angle ray modes 
were very similar to the measurements. The group delays 
measured were bracketed by curves calculated based on a 
variation of ± 10% on the solar flux about the estimated 
value for the time of measurements. 

5 . 3 Comparison of the Predictions of the Model with 

Current Techniques 

Comparisons were made of the predictions of the 
model for a path from Washington, D.C. to London, England 
for 7:00 U.T. and 23:00 U.T. and the ITS Ionospheric Pre- 
dictions for the same period. The results of these inves- 
tigations are given in Table 4. The maximum percentage 
difference between the. maximum usable frequencies for the 
paths was in all cases tested within 51. It should be 
pointed out, however, that the methods use common data for 
the predictions and the errors are not independent. The 
current prediction techniques have a long history of being 
checked against actual operating conditions, and thus the 
agreement seems to validate the results reported in this 
work, The present method does, of course, predict much 
more than just the maximum usable frequency or the optimum 
working frequency. It is upon those other parameters, such 
as the ionospheric absorption, group delays on the various 
paths, and the spatial losses, that. the importance of the 
present method rests. 
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5 . 4 Comparison of Predicted Satellite to Ground Prop- 
agation with Observations 

The diurnal variations of the electron content over 
State College, Pennsylvania and Stanford, California were 
predicted and compared with the measurements from satellites. 
The day-to-day variation was predicted and compared with the 
measurements for the different local times as the satellite 
precessed. At Stanford the predictions were compared with 
data for 4 years and 3 seasons. 

Even though the general agreement was good, discrep- 
ancies of up to 40% were shown in Figures 11, 14, and 16. 
These variations were of the same order as the day-to-day 
variation in the electron content. The predictions are for 
monthly average conditions and it remains to be determined 
whether the departures are due to changes in the neutral 
atmosphere, the ionizing radiation or some other effect. 

The assumptions of the Penn State Mk 1 Model break 
down for low latitudes, but the predictions of the electron 
content at Huancayo, Peru shown in Figure 12 were in satis- 
factory agreement. The reason for this appears to be that 
the electron content can be estimated quite well provided 
that the peak electron density and the height of the peak 
are correct, and the scale height of the ionization is well 
estimated. It is probably unwise to assume that the agree- 
ment found between the current predictions and measurements 
extends to the electron density profiles in regions far 
above the peak. 
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The gradients of the electron content were also 
predicted for a few cases. The gradients caused by predict- 
able inputs, such as sunrise and sunset can be treated well 
by the models. 

It is not possible at present to treat irregularities 
with short wave lengths such as are produced by gravity waves 
which pass through the atmosphere and produce gradients which 
change with time or by winds which are not always regular and 
can change from one night to the next. 

5 . 5 The Predictions of the Path Loss 

The predictions of the absorption at 3.36 MHz for 
the path from Annapolis to State College and at 5 MHz for 
the path from Ft. Collins to White Sands Missile Range were 
made . 

The current ionospheric models seem to give too 
little attenuation, and the electron density at the lower 
heights is probably greater than given by the Mitra-Rowe 
Model. Some models based on rocket measurements which show 
substantially higher densities were also tried, and these 
showed much greater attenuations than were observed. It 
would appear that the method can serve as a useful technique 
in checking models. 

The method showed in a very clear way the percentage 
of the absorption that took place at different altitudes. 
While it has been realized for some time that great care had 
to be exercised in interpreting ionospheric absorption 
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measurements in terms of D region densities, it was very 
difficult for people making measurements to interpret their 
data and determine what they were seeing. The present 
technique allows readily interpretable results to be pre- 
sented in a very simple manner so that an experimenter can 
determine the sensitivities of his technique or results to 
the parameters of interest. 

5 . 6 Recommendations for Future Work 

Many experimental data are now available for the 
satellite electron content measurements, oblique incidence 
ionograms, and A3 absorption measurements. Wider comparison 
of the predictions by model calculations with actual data 
are needed for different seasons and locations. 

One disadvantage of this technique is the cost of 
the ray tracing, which requires a large amount of computer 
time. In this instance, the cost of ray tracing was of the 
order of a few thousand dollars, while the cost of calculat- 
ing the complete ionospheric models was only a few dollars. 
This means that when ray tracing of this type is contem- 
plated the extra expense of a rather sophisticated iono- 
spheric model is quite small, and much greater complexities 
could be included in the ionospheric model program without 
greatly increasing the total cost. The ray tracing program 
also needs to be simplified if possible to reduce the cost 
without affecting its accuracy. 
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The Penn State Mk 1 Model and Mitra-Rowe Model must 
be coupled together by matching the boundary conditions so 
that the electron density and collision frequency profiles 
can be obtained without a discontinuity from the base of the 
ionosphere to the upper F region as a function of time, day 
number, season, and solar activity. This model coupled with 
the magnetic field model can then supply the necessary 
quantities which the ray tracing program needs to calculate 
the ray path for any desired transmitter and receiver. This 
could allow technicians to run the program by simply supply- 
ing the proper input data. In this case error estimates as 
well as parameter estimates will be required. 
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APPENDIX 1 

ESTIMATION OF INITIAL AZIMUTH ANGLE AND A GREAT CIRCLE PATH 

Before the ray tracing program can be run,it.~is 
necessary to make an estimate of the initial azimuth angle 
and of the great circle path so that the ionospheric models 
can be run for the appropriate geographical regions. 

The first parameter to be. calculated , given the 
geographic latitudes and longitudes of the transmitting and 
receiving locations, is the path length, which is taken to 
be the shorter of the great-circle distances between the 
two points, and which is computed from the spherical trig- 
onometry as shown in Figure 39. 

cos d = sin A^ sin + cos A^ cos A 2 cos (0-^ - 0 2 ) 
where 

A^ = geographic latitude of transmitter 
6 ^ = geographic longitude of transmitter 
A 2 = geographic latitude of receiver 
e 2 - geographic longitude of receiver 
d ® path length in radians 

Having obtained the path length, the bearing (the 
azimuth angle from the north pole) of the transmitter along 
the great circle path can be calculated; 
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Figure 39 Geometry of the Path 
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cos <j>i = (sin X 2 - sin X^ cos d)/cos X^ sin d 
cos <t>2 ~ (sin X ^ - sin X2 cos d)/cos X2 sin d 


where 

4 >! = bearing of transmitter to receiver in radians 
<J>2 = bearing of receiver to transmitter in radians 
Using the bearing obtained, the geographic latitude 
and longitude of any point on the great circle path which is 
a certain distance away from the transmitter can also be 
calculated in the same way. 
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APPENDIX 2 

LIST OF MAIN COMPUTER PROGRAMS 


1. Penn State Mk 1 Ionospheric Model 

This program gives the complete electron density 
profile from 210 km to 1250 km as a function of latitude, 
longitude, day number and solar activity. More details are 
described in Chapter I. This program also can calculate 
the electron collision frequency with neutrals and ions 
using Itikawa's Equation. 

2. Ray Tracing Program 

This program was originally developed by Jones (1966) 
and modified a little at Stanford Radioscience Laboratory. 

The modified version was used in this work. This program 
calculates the ray path by numerical integration of six 
simultaneous differential equations as discussed in Chapter I. 

3. Magnetic Field Model 

This program calculates the elements of the earth's 
magnetic field vector as discussed in -Chapter I. 

4. D Region Model 

This model gives the electron density profile from 
60 km to 120 km for middle latitudes (for dip latitude 30° ~ 
60°) as a function of the solar zenith angle, season, and 
solar activity. This program also gives the electron 
collision frequency with neutrals by using Itikawa's equa- 
tion and Groves ' neutral atmospheric model. 



